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ABSTRACT 

A/-methylpyrrole (Py)-A/-methylimidazole (Im) poly- 
amides are small organic molecules that bind to 
DNA with sequence specificity and can be used as 
synthetic DNA-binding ligands. In this study, five 
hairpin eight-ring Py-lm polyamides 1-5 with different 
number of Im rings were synthesized, and their 
binding behaviour was investigated with surface 
plasmon resonance assay. It was found that associ- 
ation rate (k a ) of the Py-lm polyamides with their 
target DNA decreased with the number of Im in the 
Py-lm polyamides. The structures of four-ring Py-lm 
polyamides derived from density functional theory 
revealed that the dihedral angle of the Py amide 
carbonyl is 14^18°, whereas that of the Im is signifi- 
cantly smaller. As the minor groove of DNA has a 
helical structure, planar Py-lm polyamides need to 
change their conformation to fit it upon binding to 
the minor groove. The data explain that an increase 
in planarity of Py-lm polyamide induced by the in- 
corporation of Im reduces the association rate of 
Py-lm polyamides. This fundamental knowledge of 
the binding of Py-lm polyamides to DNA will facilitate 
the design of hairpin Py-lm polyamides as synthetic 
DNA-binding modules. 

INTRODUCTION 

7V-Methylpyrrole (Py)-A^-methylimidazole (Im) poly- 
amides are small organic molecules that can recognize 



specific DNA sequences in the minor groove of B-form 
DNA, according to DNA recognition rules (1,2). Py 
favours T, A and C bases, whereas Im favours a G base. 
A lone electron pair on N-3 in Im forms a hydrogen bond 
with the 2-amino hydrogen of guanine (G). Thus, antipar- 
allel pairings of Im/Py and Py/Im specify GmC and C»G, 
respectively, and antiparallel pairings of Py/Py specify 
A»T or T»A degenerately (1,2). Aliphatic (3-alanine ((3) 
can be substituted for Py. It has been used effectively 
when molecules have more than five consecutive Py or 
Im residues, by adjusting the pitch between amide bonds 
of Py-lm polyamides and the accepting residue of the 
minor groove. Antiparallel pairings of Py/(3 and (3/Py 
specify AmT or T»A degenerately, and antiparallel 
pairings of Im/(3 and (3/Im specify G»C and C»G, respect- 
ively (3,4). 

As Py-lm polyamides can bind to DNA with sequence 
specificity comparable with DNA-binding proteins, they 
can be substituted for the DNA-binding domain of a tran- 
scription factor. Py-lm polyamides that can bind to a 
promoter region have been designed to inhibit gene ex- 
pression (5-9). Furthermore, Py-lm polyamides have 
been conjugated with a peptide or a small organic 
molecule to create synthetic transcriptional activators 
that stimulate gene expression (10-13). The dissociation 
equilibrium constants (K^s) of these Py-lm polyamides 
with their target DNA sequences were extensively 
determined, by DNase I footprinting, by the Dervan 
group. However, only a few of their corresponding asso- 
ciation rate constants and dissociation rate constants have 
been reported (14-16). It may be crucial for the design of a 
synthetic DNA-binding module to determine not only Kj^s 
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but also the association rate constant (k a ) and the dissoci- 
ation rate constant (k d ) of Py-Im polyamides because the 
association/dissociation rate constants are contingent on 
respective transcription factors (17,18). 

A zinc finger is a representative transcription factor 
DNA-binding domain that binds to GC-rich sequences 
(19). Many eukaryotic genes contain highly GC-rich se- 
quences in the promoter region (20); thus, it is potentially 
valuable to synthesize Py-Im polyamides that recognize 
the GC-rich promoter region. However, 5 -GCGC-3' and 
S'-CGCG-S' have been identified as difficult recognition 
sequences of Py-Im polyamides (3,21). It is not known 
why these sequences are difficult to recognize. 

In this study, we synthesized five 8-ring Py-Im poly- 
amides 1-5 with different numbers of Im rings, and we 
measured the k a and k d values of the DNA-Py-Im poly- 
amide complexes by surface plasmon resonance (SPR) 
assay. We also estimated the structures of four-ring Py-Im 
polyamides by density functional theory (DFT) and ab initio 
quantum chemical calculation. Our SPR data and the 
calculated structures elucidated the relationship between 
the structure of the Py-Im polyamides and the association 
rate of the Py-Im polyamides with their target DNA. 



MATERIALS AND METHODS 

General methods 

The following abbreviations apply: Fmoc, fluorenylme- 
thoxycarbonyl; DMSO, dimethyl sulphoxide; TFA, 
trifluoroacetic acid; (3, (3-alanine; y, y-aminobutyric acid; 
Py, Af-methylpyrrole; Im, Af-methylimidazole; Dp, N, 
TV-dimethyl- 1 , 3 -propylamine . 

Electrospray ionization time-of-flight mass spectro- 
metry (ESI-TOFMS) was carried out on a BioTOF II 
(Bruker Daltonics) mass spectrometer to determine the 
molecular weight of Py-Im polyamides 1-5. 

Polyamide synthesis 

Py-Im polyamides 1-5 were synthesized in a stepwise 
reaction using a previously described Fmoc solid-phase 
protocol (21). Syntheses were performed using a pioneer 
peptide synthesizer (PSSM-8, Shimadzu) with a 
computer-assisted operation system on a 36 uM scale 
(lOOmg of Fmoc- (3-alanine Wang resin). After the synthe- 
sis, Dp was mixed with the resin for 4h at 55°C and the 
mixture was shaken at 550 r.p.m. to detach the Py— Im 
polyamides from the resin. Purification of Py— Im poly- 
amides 1-5 was performed using a high-performance 
liquid chromatography (HPLC) PU-2080 Plus series 
system (JASCO), using a 10 x 150 mm ChemcoPak 
Chemcobond 5-ODS-H reverse-phase column in 0.1% 
TFA in water, with acetonitrile as eluent, at a flow rate 
of 3 ml/min and a linear gradient elution of 20-60% aceto- 
nitrile >20min, with detection at 254 nm. Collected frac- 
tions were analysed by ESI-TOFMS. 

SPR assay 

All SPR experiments were performed on a BIACORE X 
instrument at 25° C as described previously (21,22). 



The sequences of biotinylated hairpin DNAs containing 
target sequences are shown in Figure 2 and Supple- 
mentary Figure SI. The hairpin DNAs were immobilized 
on a streptavidin-coated SA sensor chip at a flow rate of 
20 jil/min to obtain the required immobilization level (up 
to ~1400 resonance units (RU) rise). Experiments were 
carried out using HBS-EP (10 mM 4-(2-hydroxyethyl)-l- 
piperazineethanesulfonic acid (HEPES), 150mM NaCl, 
3mM ethylenediaminetetraacetic acid and 0.005% 
Surfactant P20) buffer with 0.1% DMSO at 25°C, pH 
7.4. A series of sample solutions were prepared in 
HBS-EP buffer with 0.1% DMSO and were injected at a 
flow rate of 20 jil/min. To measure association and dissoci- 
ation rate constants (X D , £ a and k d ), data processing was 
performed with an appropriate fitting model using the 
BIAevaluation 4.1 program. The sensorgrams of all data 
were fitted by using the 1 : 1 binding model with mass trans- 
fer. The values of Kd, k a and k d for all data are 
summarized in Table 1 and Supplementary Table S 1 . 

Structural model calculation 

DFT and Hartree-Fock calculations were performed 
with the Gaussian 09 software package. Structural 
energy was first minimized by means of a Parameterized 
Model number 3 (PM3) model, followed by DFT and a 
6-31 1+G* polarization basis set or the Hartree-Fock 
method and a 6-3 1G* polarization basis set. 

RESULTS AND DISCUSSION 

Hairpin eight-ring Py-Im polyamide synthesis 

To investigate the binding properties of Py and Im in 
hairpin Py-Im polyamides, we designed and synthesized 
five hairpin eight-ring Py-Im polyamides 1-5 (Figure 1) by 
the Fmoc-chemistry solid-phase synthesis method. Two 
(3-alanines were attached to the N-terminal of 1-5 
(Figure 1) for the optional construction of a fluorescence 
Py-Im polyamide conjugate (Han et al., unpublished 
data). We purified 1-5 by reverse-phase HPLC, and then 
confirmed that the purity of 1-5 was >95% by analytical 
HPLC and ESI-TOFMS. The (3-Dp linker at the 

Table 1. Binding affinity of 1-5 

Py-Im polyamide a K u K a K d AG 0b 



(1(T 9 M) (10 4 M -1 s-) (1(T 2 s" 1 ) (kcal/M) 

1 ^OOOO 2.5 ± 0.50 150 ± 29 0.39 ± 0.12 -11.7 

2 (qqq2 3.7 ±2.3 120 ± 17 0.43 ± 0.25 -11.5 

3 (qqqJ 36 ± 10 28 ± 24 0 95 ± 0 70 " 10 - 2 

4 (q2o# 48 ± 4.7 30 ±4.5 1.4 ± 0.08 -10.0 

5 (£252 54 ±9.9 13 ± 2.9 0.64 ± 0.13 -9.9 



a Closed circle and open circle indicate Im and Py, respectively. 
b AG° values are calculated from the equation, AG 0 = — RTln (I/Kd), 
where R is universal gas constant, 1.987 cal/M«K; T is absolute tem- 
perature in a unit of K, here 298.15 K. 
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1 : ppPyPyPyPyyPyPyPyPypDP 




3: (3(3lmPyPyPyYPyPylmPy(3DP 




5: ppimPylmPyYlmPylmPypDP 



Figure 1. Structures of Py-Im polyamides 1-5. 



C- terminal has ~ 100-fold steric preference for A»T or T»A 
relative to GmC or CmG (23). Based on the recognition rule 
of polyamides, the target DNA sequences of 1-5 are 
5'-WWWWWW-3', 5 -WWWWCW-3\ 5 -WWWGCW-3\ 
5'-WWCGCW-3' and 5 / -WGCGCW-3 / , respectively, and 
we prepared five S^biotinylated hairpin DNAs (ODN1-5) 
(Figure 2). However, because of two P-alanines attached to 
the N-terminal of 1-5, steric hinderance between the tails of 
1-5 and the DNA minor groove may suppress the steric 
preference for A»T or TmA relative to G»C or C»G. To 
characterize the effect of the two P-alanines, we also 
prepared five 5 ; -bio tiny lated hairpin DNAs (ODN6-10) 
(Supplementary Figure SI). 

SPR assay 

To measure K u , k a and k d values of 1-5 for their target 
DNA, we performed SPR assays as described in the 
'Materials and Methods' section. The five S'-biotinylated 




2: ppimPyPyPyyPyPyPyPypDP 




4: ppimPylmPyYPyPylmPypDP 



hairpin DNAs (ODN1-5) (Figure 2) that include the target 
DNA sequences were immobilized on a streptavidin- 
coated sensor chip, and the Py-Im polyamide solutions 
were injected. As shown in Figure 3, SPR sensorgrams 
were obtained, and the kinetic binding parameters K u , 
k a and k d were determined (Table 1). The K u s of 1-5 
increased with an increase in the number of Im as 
follows: 2.5 x 10" 9 , 3.7 x 10" 9 , 3.6 x 10" 8 , 4.8 x 10" 8 and 
5.4 xlO -8 M, respectively. Interestingly, the binding 
affinities of 1-4 were 22, 15, 1.5 and 1.1-fold, respectively, 
over that of 5. Of importance was that the k d values of 1-5 
were comparable with each other (0.0039^0.014 s _1 ) 
(Table 1), whereas the k a values of 1-5 were 1.5 x 10 6 , 
1.2 x 10 6 , 2.8 x 10 5 , 3.0 x 10 5 and 1.3 x 10 5 M -1 s"\ re- 
spectively (Table 1). These results indicate that the asso- 
ciation rate of the Py-Im polyamides with their target 
DNA decreased as the number of Im in the Py-Im poly- 
amides increased. However, once 1-5 bound to their target 



Nucleic Acids Research, 2012, Vol. 40, No. 22 11513 



0DN1 5 3, 



0DN2 



0DN3 



0DN4 



0DN5 



-XGCCATATATGGC T 
-CGGTATATACCG T 

T 



-XGCCATATCTGGC T 
- C GGT ATAGAC C G _ T 
T 



-XGCCATAGCTGGC T 
-CGGTATCGACCGT 

T 



-XGCCATCGCTGGC T 
-CGGTAGCGACCG T 
T 



-XGCCAGCGCTGGC T 
-CGGTCGCGACCGT 

T 



Figure 2. Sequences of S'-biotinylated hairpin DNAs (ODN1-5). 
X represents biotin. The binding sequences of the Py-Im polyamides 
are shown in bold. 



DNAs, the dissociation rates of the Py-Im polyamides 
from the respective complexes were comparable with 
each other. We also determined the free energy change 
(AG°, kcal/M) from the on the formation of the 
Py-Im polyamides 1-5/DNA complexes (Table 1). The 
AG° values of 1/ODN1, 2/ODN2, 3/ODN3, 4/ODN4 
and 5/ODN5 were -11.7, -11.5, -10.2, -10.0 and 
— 9.9kcal/M, respectively. 

We also prepared five S'-biotinylated hairpin DNAs 
(ODN6-10) (Supplementary Figure SI) to characterize 
the effect of the two (3-alanines, and as shown in Supple- 
mentary Figure S2 and Supplementary Table SI, we 
determined the kinetic binding parameters k a and 
k d . As described earlier in the text, the ^ D s of 1-5 also 
increased with an increase in the number of Im as follows: 
6.1xl0~ 9 , 1.3xl0~ 8 , 2.5xl0~ 8 , 1.1 x 10~ 7 and 
2.5 x 10 -7 M, respectively, and the binding affinities of 
1-4 were 41, 19, 10 and 2.3-fold, respectively, over that 
of 5. The binding affinities of 1-5 to ODN6-10 were 
reduced 2.3, 3.5, 0.72, 2.3 and 4.6-fold, respectively, 
compared with those to ODN1-5 (Table 1 and 
Supplementary Table SI), indicating that the [3-Dp 
linker at the C-terminal, with two (3-alanines attached to 
the N-terminal, has a slight steric preference for A»T or 
T»A relative to G»C or C»G. The k d values of 1-5 were 
also comparable with each other (0.01 1~0. 023 s _1 ) 
(Supplementary Table SI), whereas the k a values of 1-5 
were 1.8 x 10 6 , 1.8 x 10 6 , 5.1 x 10 5 , 1.8 x 10 5 and 8.9 x 10 4 
M _1 s _1 , respectively (Supplementary Table SI). The AG° 
values of 1/ODN6, 2/ODN7, 3/ODN8, 4/ODN9 and 
5/ODN10 were -11.2, -10.8, -10.4, -9.5 and 
—9.0 kcal/M, respectively. 

Previously, Crothers and coworkers reported the £ a and 
k d values of ImPyPyPy-y-ImPyPyPy-P-Dp, and the k d 
value was 0.002 s _1 , which is consistent with our data. 
However, the k a value was 7.0 x 10 7 M _1 s _1 and 46-fold 
times higher than that of 1. Dervan and coworkers (24) 
suggested that the Im located at the C-terminal end of 
each four-ring Py-Im polyamide subunit is somehow less 
capable of strong hydrogen bond formation than the 
N-terminal residues. Therefore, the binding affinity of 



ImPyPyPy-y-ImPyPyPy-P-Dp may be relatively high, 
like that of ImlmPyPy-y-ImlmPyPy-P-Dp as discussed 
later in the text. 

Structures of four-ring Py-Im polyamide estimated by 
density functional theory 

To obtain insight into the differences among the observed 
Kd values, we calculated the model structures of five 
4-ring Py-Im polyamides, PyPyPyPy, ImPyPyPy, 
PyPylmPy, ImPylmPy and Imlmlmlm by DFT, as 
described in the 'Material and Methods' section. Py-Im 
polyamide 1 contains two PyPyPyPy, 2 contains 
ImPyPyPy and PyPyPyPy, 3 contains ImPyPyPy and 
PyPylmPy, 4 contains ImPylmPy and PyPylmPy and 
5 contains two ImPylmPy. For comparison, the structure 
of Imlmlmlm was also calculated. The optimized struc- 
tures and selected torsion angle of four-ring Py-Im are 
shown in Figure 4. 

It was found that PyPyPyPy has a helical structure, and 
the dihedral angles of the Py carbonyl of the amide 
N1-C2-CO-0 are ~17°, with a slightly smaller dihedral 
angle at the C-terminal Py. This large dihedral angle of the 
Py residue is derived from a large steric hindrance between 
the H of C-3 in Py and the H of the contiguous amide. In 
clear contrast, the corresponding dihedral angles of the Im 
residue in the four-ring Py-Im tetramers are small (<l-2°) 
because of the lack of steric hindrance, and the planarity 
of the molecules is increased. As a result, Imlmlmlm 
becomes almost planar (Figure 4). Because of the double 
bond between C2 and N3 in the Im, a lone electron pair of 
N3 face a H of the contiguous amide, and a hydrogen 
bond may form between N3 in the Im and H in the con- 
tiguous amide, resulting in the shorter distance of C-N in 
the amide and between two Ns in the contiguous amides 
and the more acute angle of the two ring-to-amide bonds, 
compared with Py, and the coplanarity of the Im ring 
and the contiguous amide (Supplementary Figure S3). 
The angles of the two ring-to-amide bonds in each Py 
and Im ring were ~147° and 138°, respectively 
(Supplementary Figure S3), which is consistent with a 
previous result (25). 

We also determined the N-to-N distances at both ends 
of each of the four-ring Py-Im polyamides. The N-to-N 
distances for PyPyPyPy, ImPyPyPy, PyPylmPy, 
ImPylmPy and 0 lmlmlmlm were 15.59, 14.83, 13.72, 
12.91 and 10.64 A, respectively (Figure 4). We also calcu- 
lated the four-ring structures by ab initio quantum 
chemical calculation. The calculated structures were 
almost the same as those derived by DFT (Supplementary 
Figure S4). The calculated structural data suggest that 1 is 
the least curved Py-Im polyamide, and 5 is the most 
curved Py-Im polyamide among 1-5. 

Implication of binding of hairpin Py-Im polyamide to 
DNA minor groove 

It has been pointed out by the Dervan group that the 
twisted shape of Py-Im polyamides fits well into the 
minor groove of DNA. In a DNA-cyclic polyamide 
complex, relatively large torsion angles were observed in 
Im residues and Py residues (26,27). These results clearly 
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Figure 3. SPR sensorgrams for the interaction of Py-Im polyamides with hairpin DNAs immobilized on a sensor chip SA. (A) Py-Im polyamide 1 
with ODN1 at a concentration range from 1.25 x 10~ 8 (lowest curve) to 2.0 x 10~ 7 M (highest curve). (B) Py-Im polyamide 2 with ODN2 at a 
concentration range from 1.25 x 10~ 8 (lowest curve) to 2.0 x 10~ 7 M (highest curve). (C) Py-Im polyamide 3 with ODN3 at a concentration range 
from 2.5 x 10~ 8 (lowest curve) to 4.0 x 10~ 7 M (highest curve). (D) Py-Im polyamide 4 with ODN4 at a concentration range from 2.5 x 10~ 8 (lowest 



curve) to 8.0 x 10" 
(highest curve). 



M (highest curve). (E) Py-Im polyamide 5 with ODN5 at a concentration range from 2.5 x 10 (lowest curve) to 8.0 x 10 M 



indicate that a large conformational change is necessary 
for Im residues on DNA binding. The fact that the angles 
of the two ring-to-amide bonds in each Py and Im ring 
were ~147° and 138°, respectively, which indicates that Py 
has less curvature, and Im has too much curvature. This 
also suggests that the conformational change of Im 
residues could lead to a match to the regular B-form 
DNA, and 5 requires more energy than the other four 
Py-Im polyamide 1-4 to change the structure for 
binding to the target DNA, resulting in the slowest asso- 
ciation rate of 5 with the target DNA among 1-5. 

Dervan and coworkers also reported other Py-Im poly- 
amides that recognize the GC-rich sequences 5 / -GGGG-3 / 
and 5-GGCC-3' (24,28). In the case of the Py-Im poly- 
amide that recognizes 5 -GGGG-3', the association equi- 
librium constant (K A ) is 2.8 x 10 7 IVT 1 (The K u is 



3.6 x 10~ 8 M) (24). The K A of Py-Im polyamides that 
recognize 5 / -GGGW-3 / , 5 / -GGWW-3 / and 5 / -GWWW-3 / 
have been determined as 3.7 x 10 8 , 5.0 x 10 8 and 3.5 x 10 9 
M _1 , respectively (28-31). Among these four Py-Im poly- 
amides, the association equilibrium constant decreased as 
the number of Im in the Py-Im polyamides increased. 
However, the K A of the Py-Im polyamide that recognizes 
S'-GGCC^' was 9.7 x 10 9 IVT 1 (24), which is higher than 
that of the Py-Im polyamide that recognizes 
y-GWWW-3'. The K A of the Py-Im polyamide that rec- 
ognizes S'-GGCW^ was determined to be 4.0 x 10 8 M -1 
(32), which is similar to that of 5 / -GGGW-3 / . The Py-Im 
polyamide that recognizes 5 -GGCC-3' may be an excep- 
tional Py-Im polyamide, like ImPyPyPy-y-ImPyPyPy-P- 
Dp that recognizes 5 / -GWWC-3 / as described previously. 
However, the association equilibrium constant of the 
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1.20° 

Figure 4. Models of four-ring polyamide subunits calculated by density functional theory. Angles indicate the dihedral angles between each Py or Im 
and the contiguous amide (N1-C2-CO-0). (A) Py-Py-Py-Py. Dotted line indicates the N-to-N distance at both ends of the four-ring polyamide. 
(B) Im-Py-Py-Py, (C) Py-Py-Im-Py, (D) Im-Py-Im-Py and (E) Im-Im-Im-Im. 



other three Py-Im polyamides decreased as the number of 
Im in the Py-Im polyamides increased. Our calculated 
structural data indicate that ImPyPyPy is less curved 
compared with PyPylmPy (Figure 4B and C). These 
results suggest that not only the number of Im, but also 
the position of Im influences the Py-Im polyamide struc- 
ture and the association rate of the hairpin eight-ring 
Py-Im polyamide to the target DNA. 



As reported previously, replacement of Py with an ali- 
phatic (3-alanine can increase binding affinity and provide 
flexibility in the polyamide structures, and the binding 
affinity of Im-(3-ImPy-y-Im-(3-ImPy-(3-Dp that recognizes 
5 / -GCGC-3 / was 100-fold over that of ImPylmPy-y- 
ImPyImPy-(3-Dp (3). Measurement of Kj), k a and k& 
values of Py-Im polyamides containing Py/P and/or 
Im/P pairs is important for the next step of Py-Im 
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polyamide design. We have also replaced two Py in 5 by 
(3-alanine, resulting in construction of (3-(3-Im-(3-ImPy-y- 
Im-(3-ImPy-(3-Dp, and we measured the Kd, k a and k d of 
(3-(3-Im-(3-ImPy-y-Im-(3-ImPy-(3-Dp. Interestingly, the k a 
and k d of (3-(3-Im-(3-ImPy-y-Im-(3-ImPy-(3-Dp were 
improved by ~ 10-fold, compared with those of 5 for 
ODN5 or ODN10 (Y.-W. Han et al., unpublished data). 
Further analysis of Py-Im polyamides containing Py/P 
and/or Im/(3 pairs is now in progress. 

CONCLUSION 

In this study, using SPR assays, we measured the k a and 
k d of Py-Im polyamides 1-5 to characterize Py and Im in 
hairpin Py-Im polyamides in more detail. Because k a and k d 
of some transcription factors have been determined and 
were contingent on the respective transcriptional factors, 
the measurement of k a and k d of Py-Im polyamides is 
also crucial for the design of a Py-Im polyamide as a syn- 
thetic DNA-binding module of a transcription factor. SPR 
data demonstrated that the k d values of 1-5 were between 
0.0039 and 0.014 s" 1 . The k a values of the Py-Im poly- 
amides decreased as the number of Im in the Py-Im poly- 
amides increased. DFT calculations suggest that an increase 
in planarity, induced by the incorporation of Im, reduced 
the association rate of Py-Im polyamides. These data 
indicate that the number and the position of Im in Py-Im 
polyamides influence the k a but not the k d of the Py-Im 
polyamides; thus, enabling us to estimate the DNA 
binding kinetics of Py-Im polyamides. 

We synthesized Py-Im polyamide 1-5 which contained 
two (3-alanines at the N-terminal in this study, and SPR 
data also demonstrated that the (3-Dp linker at the 
C-terminal of 1-5 had a slight steric preference for A»T 
or T»A relative to G»C or C»G. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online: 
Supplementary Table 1 and Supplementary Figures 1-4. 
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